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Abstract 

A generalization of the notion of ellipsoids to curved Riemannian spaces 
is given and the possibility to use it in describing the shapes of rotating bod- 
ies in general relativity is examined. As an illustrative example, stationary, 
axisymmetric perfect-fluid spacetimes with a so-called confocal inside ellip- 
soidal symmetry are investigated in detail under the assumption that the 
4-velocity of the fluid is parallel to a time-like Killing vector field. A class 
of perfect-fluid metrics representing interior NUT-spacetimes is obtained 
along with a vacuum solution with a non-zero cosmological constant. 

PACS numbers: 04.20.Jb, 04.40.Nr 



1 Introduction 

In the past much effort has been made to find an exact solution of Einstein's field 
equations representing the gravitational field of a compact rapidly rotating object, 
but without much success. In particular, a great deal of work has been done 
concerning stationary, axially symmetric perfect-fluid configurations as possible 
candidates for modelling the interior of rapidly rotating bodies. In addition, in 
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order to give a global space-time model completely describing the gravitational 
field of the source, the possibilities for matching an interior perfect-fluid solution 
to an exterior vacuum metric have been extensively studied. So far, however, the 
only known physically acceptable global exact solution in the examined category 
is the so-called Neugebauer-Meinel disk of dust [T], which may be considered as 
the limiting case of uniformly rotating relativistic objects. 

Apart from the above mentioned disk, in many cases encountered in the search 
for physically acceptable solutions it is not possible to give an obvious geometrical 
interpretation for the discovered configurations. On the other hand, it is well 
known that in Newtonian gravity numerous results have been obtained for the case 
in which a rotating fluid in equilibrium has an ellipsoidal shape [2|. In the present 
work we will constrain our attention to those rapidly rotating configurations in 
general relativity, for which there is a natural way for imposing the restriction that 
the surface of the considered object is ellipsoidal. Using the approach presented in 
this paper the matching surface between the interior and exterior region may be 
conveniently chosen to be the ellipsoid representing the boundary of the rotating 
fluid body. It should be emphasized, however, that other shapes for a rotating 
fluid might also be possible in both general relativity and Newtonian theory as 
well. 

In the Euclidean geometry of flat space the notion of an ellipsoid is clear and 
easily understood. In curved spaces, however, it does not have such an obvious 
intuitive meaning. Therefore, in the general-relativistic analysis of perfect-fluid 
bodies having an ellipsoidal shape, it is necessary to determine first what an 
ellipsoid in curved space is. Originally, it has been proposed by Krasihski |3] to 
consider a 3-dimensional, axially symmetric Riemannian space in which the line 
element ^'ds 2 of the 3-dimensional space can be written as 

^ds 2 = z(j), 6)dp 2 + (p 2 + a 2 cos 2 6)d6 2 + (p 2 + a 2 ) sin 2 9d<p 2 , (1.1) 

where a is a constant and z is a positive function. In such spaces the coordinate 
surfaces p = const represent the analogue of a one-parameter congruence of 
confocal Euclidean ellipsoids of revolution. Later a more general, coordinate-free 
definition of (not necessarily axially symmetric) Riemannian spaces that can be 
filled in with a one-parameter family of confocal ellipsoids was given by Racz 
[3], recovering the metric (jl.lj) as a special case. In [3| it was assumed that the 
union of local rest spaces of the observers co-moving with matter has the metric 
([1.1)1 . whereas in jl] it was the 3-space of timelike Killing orbits which had this 
ellipsoidal symmetry. By using the approach of 4 , a uniqueness theorem for 
stationary axisymmetric vacuum spacetimes with a particular type of ellipsoidal 
symmetry was proven in [5]. We remark here that ellipsoidal shapes in curved 
spaces, but in a different context, have been considered in || as well (see also 
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references therein). 

The purpose of the present paper is to construct a framework for considering 
ellipsoidal shapes in general relativity, which covers the situations studied in [3] 
and [1] as special cases. In this way a firm mathematical basis for the ideas 
pioneered by Krasihski [3 is given and the constructions introduced by Racz jl] 
are extended in order to include a wider class of spaces in the analysis. This 
work also aims to examine the applicability of the constructed framework on the 
particular example of stationary, axially symmetric rigidly rotating perfect fluid 
spacetimes with certain ellipsoidal symmetries as test models for the interior 
solution. It has been shown by Geroch and Lindblom [7j that rigidly rotating 
perfect fluids can be considered to be the equilibrium configurations of relativistic 
dissipative fluids so they are of obvious astrophysical interest. Although the 
test models considered in this paper did not provide perfect-fluid solutions with 
physically suitable properties, the existence of such solutions within the very wide 
class of spacetimes with some kind of ellipsoidal symmetry cannot be excluded. 

This paper is organized as follows. First the definition of a Riemannian 3-space 
that can be filled in with a one-parameter congruence of concentric ellipsoids is 
given. In such curved spaces one can introduce a notion analogous to the notion 
of ellipsoids in Euclidean space. Next, the possibilities for considering this type of 
3-spaces within the framework of general relativity are discussed. The definition 
of a "conformally ellipsoidal spacetime" as well as the more restrictive defini- 
tions of a "spacetime with co-moving ellipsoidal symmetry" and of a "spacetime 
with inside ellipsoidal symmetry" are presented. As an illustration of these no- 
tions, stationary, axially symmetric, rigidly-rotating perfect-fluid configurations 
with "confocal" inside ellipsoidal symmetry are considered. After reviewing some 
basic notions pertaining to rigidly rotating fluids, our considerations are first spe- 
cialized to the subcase in which the gradient of the norm of the time-like Killing 
vector field and its twist are linearly dependent. Making use of the integrability 
conditions of the field equations, the general form of the metric for the selected 
perfect-fluid configurations is determined and it is demonstrated that the result- 
ing metric can be joined smoothly to the Taub-NUT solution. Some previously 
known exact solutions are shown to belong to the class of spacetimes studied 
here. In the complementary case when the norm of the gradient and the twist of 
the Killing field are linearly independent it is shown that a particular (possibly 
new) perfect-fluid solution satisfies the integrability conditions. The resulting 
particular metric, however, can be re-interpreted as a vacuum solution with a 
non-zero cosmological constant. Finally, the Appendix contains a review of some 
facts concerning ellipsoids in flat Euclidean space. 
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2 "Ellipsoids" in curved space 

In the following a possible approach for considering ellipsoidal shapes in curved 
Riemannian spaces is presented. In particular, the definition given in jl] is gen- 
eralized and it is shown how a convenient coordinate system can be chosen in the 
Riemannian spaces under consideration. 

Definition 1: A Riemannian 3-space (<S, h^) can be filled in with a one-parameter 
congruence of concentric ellipsoids if it satisfies the following properties: (a) there 
exists in it a congruence of surfaces, £ p , with induced metric Xij[p]i (b) there exists 
a diffeomorphism, mapping S into a subset, ^[S], of R 3 so that the tensor 
field is just the induced metric field for some one-parameter congruence 

of concentric ellipsoids in ^[S] C R 3 and (c) for any value of p, the extrinsic 
curvature tensor field, /%[/?], of the surfaces £ p , embedded in S, is conformal to 
the metric tensor of the Euclidean unit sphere, 7^. 

We will refer to the surfaces £ p as "ellipsoids in curved space" . Note that this 
definition has been obtained by generalizing the properties of Euclidean ellipsoids 
(given in the Appendix) to curved spaces, just as it has been done in However, 
an essential difference between the above definition and the one given in 0] is 
that the term "confocal" used in 0j is replaced here by the term "concentric". 
This allows one to consider Riemannian spaces of a much more general character 
than it was the case in 

By analogy with the 3+1 decomposition of a spacetime (A4,g a b) (see, e.g., 
§21.4 of jH]), a 3-dimensional Riemannian space (S, h a b) may be decomposed in a 
similar way in 2+1 dimensions. Let £ p denote a surface labelled by the parameter 
value "p" in a one-parameter congruence of 2-dimensional surfaces embedded in 
a 3-dimensional Riemannian space. Since there is always only one surface £ p of 
the congruence passing through any given point of S, the parameter p used for 
labelling the surfaces may be chosen as one of the coordinates of S. In such 
a coordinate system the distance between two nearby points belonging to two 
different surfaces £ p and £ p +dp can be given as 



where a and (3 l are the 2+1 analogues of the "lapse function" and "shift vector", 
Xij[p] is the induced metric on £ p , while x l = 2, 3) is a coordinate system on 
this surface. Expression ()2.1j) gives, in fact, the line element of the 3-dimensional 
Riemannian space (S,h a b)- The components of the metric hab and its inverse, 
h ab , are given by 



ds 2 = a 2 dp 2 + x%j\p]{dx i + (3 t dp)(dx j + (3 j dp), 
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1/a 2 -ft/a 2 
-P i /a 2 x ij + Pft/oi 2 



(2.3) 



The requirement that ^/*Xij [p] is the induced metric field for some one-para- 
meter congruence of concentric ellipsoids in ^[S] C R 3 implies that there exists 
a coordinate system (x 2 = 9, x 3 = 0) in which the line element of the induced 
metric on the surface S p can be written as (compare with equation (|A.2J) of the 
Appendix) 



ds2 {Xij[p\) 



K x cos 



+ Uq sin 2 0) cos 2 9 + p 2 sin 2 9 



d9 2 



+ (x Q sin <p + y cos <p) sin 9 ■ d(f) , 
+ 2(i/q — Xq) sin 9 cos 9 sin </> cos (p ■ dOdcf) 



(2.4) 



where xq = xo(p) and yo = yo(p). On the other hand, the requirement that the 
extrinsic curvature tensor field Kij[p] is conformal to the metric of the Euclidean 
unit sphere, 7^, means that there exists a smooth (at least C 2 ) everywhere posi- 
tive function k(p, 9, (p) so that 1 



1 v 



Hp, 4>)n. 



(2.5) 



In general, the extrinsic curvature, or "second fundamental form", of a surface 
is given by the expression [HI [TU] 



dx a dx b 
WW 



D a n b 



(2.6) 



where n a is the unit normal to the surface, D a is the covariant derivative as- 
sociated with the metric of the Riemannian 3-space, x a are coordinates on the 
3-space, while are coordinates on the 2-dimensional surface (a, b = 1,2,3 and 
i,j = 2, 3). If the surface is given by an equation of the form 



then the unit normal can be given as 







± 



h 



bc df Of 



dx b dx c 



-V 2 df 
dx a 



(2.7) 



{21 



with the ambiguity in sign arising from the ambiguity in the orientation of the 



normal. In the present case x l 



p, x' 



9 or 



the surface 



1 From now on we omit the labels "[/>]" in the notation and denote any tensor field Tij[p] 
simply as r%y Thus, the dependence on p, as well as on the coordinates 9 and <j), is understood. 
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equation is given by / = p — const = 0, and, thus, the components of the unit 
normal are (ni,n 2 ,n 3 ) = (±a,0,0). 

In this setup, with the aid of (|2.2|) and ()2.3|) one obtains from ()2.6|) 

for the extrinsic curvature (see also [Hj), where T^- are the Christoffel symbols for 
the metric Xiji which is given by (|2.4j) . Furthermore, in the coordinate system of 
(12.4)1 . the line element of the unit sphere may be written as 

ds 2 ^) = d9 2 + sin 2 9dtf, (2.10) 

so it follows from equation ()2.5|) that 

= sin 2 fl, « w) = 0, (2.11) 

where K(e<f>) and K(^) each denote the corresponding component of the 

extrinsic curvature. Since for the spaces under consideration the metric Xij has 
been already specified by equation (|2.4jl . it follows from (J2.9)) that the system 
(|2.11j) contains as unknown functions only the two components of 

Now it is straightforward to show that (|2.11j) is identically satisfied if one 
choses the components (3^) and (3^) of $ as 

= [(^oXq cos 2 + yoy'o sin 2 — p) cos9 + F(p)] sin9, (2-12) 
= [(— xox' + yoy'o) sin cos + G(p)]sin 2 9, (2-13) 

where F(p) and G(p) are arbitrary functions and the prime denotes derivation 
with respect to p. Note that (3^) and f3^ satisfy the system ()2.11|) for arbitrary 
choice of the functions F and G, including the particular case F = G = 0. In fact, 
fixing the functions F and G merely corresponds to a definite choice of coordinate 
system. Therefore, without loss of generality, hereafter these two functions will 
be set to be identically zero. Correspondingly, there exists a coordinate system 
in which the line element of a curved Riemannian 3-space that can be filled in 
with a one-parameter congruence of concentric ellipsoids takes the form 

ds 2 = z (p, 9, (f)) ■ dp 2 + [(xl cos 2 (j) + y% sin 2 4>) cos 2 9 + p 2 sin 2 9} ■ d9 2 

+ (xq sin 2 + y 2 cos 2 0) sin 2 9 ■ d(f> 2 

+ 2(x x' cos 2 + y o y' sin 2 — p) cos^sin^ • dpd9 , (2-14) 

+ 2(— x x' + yoy'o) cos sin sin 2 9 ■ dpd<p 

+ 2(?/q — Xq) cos sin cos 9 sin 9 ■ d9d<p 
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where z(p, 9, 0) is a C 2 function, x and y are functions depending merely on the 
coordinate i p\ and the prime denotes derivation. In these coordinates the p = 
const surfaces are the curved-space analogues of Euclidean ellipsoids. Note also 
that the particular choice z = ((x' ) 2 cos 2 <t>+{y'o) 2 sin 2 0) sin 2 #+cos 2 9 recovers the 
metric of flat Euclidean space in ellipsoidal coordinates (compare with equation 
(lA~2l) of the Appendix). 

The freedom carried by the functions x and y as well as the ^-dependence 
of z in the above form of the metric implies that (|2.14|) describes a much wider 
class of spaces than the metric (J1.1J) . Putting, e.g., Xq = p 2 + a 2 , y 2 = p 2 + b 2 into 
(|2.14|) where a and b are constants, recovers the line element of non-symmetric 
Riemannian spaces that can be filled in with a congruence of confocal ellipsoids, 
originally given in [I] 2 . The present paper, however, will be particularly concerned 
with the special case in which z(p, 9, (p) = z(p, 9) and %l = y 2 = p 2 + a 2 , where 
a = const. Then the metric of the considered 3-dimensional Riemannian space 
reduces exactly to the metric (jl.lj) . In this special case the congruence of the p = 
const surfaces is the analogue of a congruence of rotationally invariant confocal 
oblate ellipsoids in Euclidean space. 

3 Ellipsoidal shapes in general relativity 

In the previous section we have been concerned merely with structures in 3- 
dimensional Riemannian geometry, without connecting them to any structures 
in general relativity. In this section Riemannian spaces within the 4-dimensional 
spacetime of general relativity will be considered, assuming that these Rieman- 
nian spaces can be filled in with a one-parameter family of concentric ellipsoids, 
in the sense of Definition 1 from the previous section. 

In general relativity a description of any shape depends on the motion of the 
observers performing the description, as it has been already emphasized in [3]. 
E.g., a body that has an ellipsoidal shape as seen by one class of observers might 
appear to have a completely different shape as seen by some other observers. 
Therefore when considering ellipsoidal shapes in general relativity, one should 
also specify a choice of preferred observers. To start off, consider the following 
construction: 

Definition 2: A spacetime {M,g a ^) is conformally ellipsoidal if there exists in 
it a congruence of preferred observers for which the union of their local rest 

2 Note, however, that the non-diagonal term with "dOdcj)" is missing form the line element 
given by equation (11) in py. Nevertheless, this omission does not influence the remaining 
results in |]y, because they correspond to the rotationally invariant case, in which the mentioned 
non-diagonal term identically vanishes. 
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spaces is conformal to a Riemannian 3-space (<S, h ab ) that can be filled in with a 
one-parameter congruence of concentric ellipsoids. 

For the observers moving with 4- velocity u a , Definition 2 implies 



where Q 2 is a conformal factor. A spacetime is ellipsoidal if Q 2 = 1. Note that 
the original term "ellipsoidal spacetime" first introduced in j3] is equivalent to 
the notion of a stationary, axisymmetric spacetime with "co-moving ellipsoidal 
symmetry" (see Definition 3) in the terminology of the present work. 

Notice now that the choice of preferred observers has still not been uniquely 
specified. Therefore, in accordance with the remark above Definition 2, in order 
to consider ellipsoidal shapes within general relativity, Definition 2 should be 
supplemented by a prescription for choosing the preferred observers. There is an 
enormous freedom in fixing the congruence of observers and the conformal factor 
referred to in Definition 2. In particular, we will consider two possible situations, 
which give rise to the definitions presented bellow, generalizing the constructions 
given in |3] and 0j. 

Let us recall that the main motivation for considering the analogues of ellip- 
soidal surfaces in curved spaces was examining rapidly rotating bodies having an 
ellipsoidal surface in the framework of general relativity. As discussed in |3] , a nat- 
ural 3-dimensional Riemannian space in which one should consider the ellipsoidal 
shapes of rotating bodies in general relativity is the union of local rest spaces of 
observers co-moving with matter. This serves as a motivation for Definition 3. 

Definition 3: A non-empty spacetime possesses co-moving ellipsoidal symmetry 
if it is ellipsoidal (i.e., conformally ellipsoidal with Q 2 = 1) and if the class 
of preferred observers is exactly the class of observers co-moving with matter. 
Furthermore, any ellipsoidal vacuum spacetime will be said to possess co-moving 
ellipsoidal symmetry, by definition. 

In an empty spacetime, in general, there is no way to specify a class of distin- 
guished observers. However, if there is a Killing vector field present on spacetime, 
one may give, for instance, the following definition, regardless whether the space- 
time is empty or non-empty: 

Definition 4: A spacetime with a timelike Killing vector field £ a possesses in- 
side ellipsoidal symmetry if it is conformally ellipsoidal with Q 2 = (— £ a £ a ) _1 =: 
(— v ) _1 and if the 4-velocity of the preferred observers, u a , is a unit vector field 
aligned with the Killing vector field, i.e., if u a = (— t>)~ 1//2 £ a . 

In particular, if the rest spaces of the preferred observers from Definition 4 can 
be filled in with a congruence of confocal ellipsoids, we will refer to the spacetime 



9ab + u a u b = Q h ab , 
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as possessing confocal inside ellipsoidal symmetry* . 

A spacetime with co-moving ellipsoidal symmetry need not, but might possess 
inside ellipsoidal symmetry and vice versa. For instance, the Kerr spacetime is 
an example of a spacetime which possesses both these types of ellipsoidal sym- 
metry at the same time. The present paper will be mainly concerned with the 
case of inside ellipsoidal symmetry, while spacetimes with co-moving ellipsoidal 
symmetry are the subject of ongoing research. It turns out, however, that a class 
of metrics with inside ellipsoidal symmetry discussed in this paper possesses co- 
moving ellipsoidal symmetry as well (see equations ()4.20j) and (j4.21|) later on in 
this paper). 

In jH] stationary, axially symmetric empty spacetimes with confocal inside 
ellipsoidal symmetry have been considered 4 . By examining the integrability con- 
ditions of the field equations it has been demonstrated that the full set of solutions 
to the vacuum Einstein equations satisfying these conditions coincides with the 
family of the Kerr-NUT solutions with vanishing electric and magnetic charges. 
In the remaining part of the present paper a method similar to that of will 
be used to study stationary, axisymmetric perfect-fluid configurations possessing 
confocal inside ellipsoidal symmetry and having the 4- velocity of the fluid parallel 
to a timelike Killing vector field. 

4 Rigidly rotating perfect-fluid spacetimes with 
confocal inside ellipsoidal symmetry 

4.1 The field equations for rigidly rotating perfect fluids 

In the following we will restrict our attention to stationary, axisymmetric perfect- 
fluid spacetimes with the energy momentum tensor 

Tab = (ju + P)u a u b + Pg ab , (4.1) 

where /i and P are the density and pressure of the fluid and for which the 4- 
velocity of the fluid is a unit vector field u a aligned with a timelike Killing field 
£ a , so that 

u a = (-^r l/2 r, (4.2) 

where v = £ a £ a is the norm of the time-like Killing vector field £ a . Such fluids are 
"rigid" in the sense that they are expansion- and shear-free. Moreover, for the 

3 We remark that the term "inside ellipsoidal symmetry" originally introduced in 0j is equiv- 
alent to the term "confocal inside ellipsoidal symmetry" used in the present work. 

4 Note that [5] refers to "ellipsoidal spacetimes" which are, in fact, spacetimes with confocal 
inside ellipsoidal symmetry in the terminology of the present paper. 
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selected configurations the twist of the Killing field, defined by u a : = £ a fccd£ 6 V c £ d , 
can be given as a gradient of a function, that is, u a = V a co>, where to is called the 
twist potential. The existence of the Killing field allows one to use the projection 
formalism of general relativity (for a review see, e.g., ^UIHEIJE]) anc ^ ^° WT ^ e 
the field equations for the selected perfect-fluid configurations in the form [13 HE] 

R ab - 1671V- 1 Ph ab = V 2 [(D a v) (D b v) + {D a u) {D b u)\ , (4.3) 

D m D m v = v- l [{D m v)(D m v) - (D m u){D m u)\ - 8tt(/x + 3P) , (4.4) 

D m D m u = 2v~ 1 (D m u)(D m v) , (4.5) 

where R ab and D a are the Ricci-tensor and covariant derivative associated with 
the 3-dimensional Riemannian metric h ab on the space of time-like Killing orbits 
defined by 

{-vY l h ab = g ab f a f 6 . (4.6) 

v 

In addition to the above equations, the equation of motion 

9a P + hv + P)^ = 0, (4.7) 
2 v 

which follows from the above ones, is also satisfied. 

In the spirit of Definition 4, for spacetimes with confocal inside ellipsoidal 
symmetry the 4-velocity of the preferred observers is given by u a = (-v)~ 1/2 C, 
i.e., in the present case of rigidly rotating fluids it coincides with the 4-velocity of 
the fluid elements. Therefore, the preferred observers referred to in Definition 4 
are now co-moving with the fluid. This, however, does not necessarily mean that 
the spacetime possesses co-moving ellipsoidal symmetry in the sense of Definition 
3, because that would also require fl 2 = 1 for the conformal factor Q 2 introduced 
in equation ([3.1)1 . In the present case fl 2 = (— f) _1 and h ab is expected to be 
the metric of an axially symmetric Riemannian 3-space that can be filled in 
with a one-parameter congruence of confocal ellipsoids. Thus, equation ([3.1)1 is 
equivalent now to equation ([4.6)1 . which defines the metric on the space of timelike 
Killing orbits. Using p. 1)1 for the line element of the metric h ab , it follows that in a 
coordinate system (x° = t, x 1 = p, x 2 — 8, x 3 — 0) 5 adapted to the Killing vector 

5 Note that, in general, the coordinates x , x 2 and x 3 of spacetime may depend on the 
coordinates p, 9 and <f> of the three-dimensional Riemannian space in a more complicated way. 
In the present case of rigidly rotating fluids with inside ellipsoidal symmetry, however, the 
simple identifications x 1 — p, x 2 = 9, x 3 — <fi are possible. 
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fields, the line element of a stationary, axially symmetric spacetime possessing 
confocal inside ellipsoidal symmetry can be given in the form 



ds 2 = v{p,9){dt + A(p,9)d(j)f 

- —, — -r (z(p, 6)dp 2 + (p 2 + a 2 cos 2 6)d6 2 + (p 2 + a 2 ) sin 2 9d(jA , (4.8) 
v(p, 9)\ J 

where v(p,8), A(p,8) and z(p,8) are the unknown metric functions. The metric 
function A(p, 9) is recovered from the quantities used in the projection formalism 
through the equations 

dA Ju doj dA p 2 + a 2 doo . . 

op v 2 o9 69 v l ^u op 

where we have introduced the function 

2 2 

u = u(p, 9) := z(p, 9) / \ a (4.10) 
p z + a z cos z 9 

In the coordinate system of ()4.8|) . from the field equation ()4.3|) one has 

16™-^=— (4.11) 

and therefore one can write 

Hab ■= Rab ~ ^h AB = \v- 2 [{d A v){d B v) + (d A Lu)(d B cu)} , (4.12) 

where A, B — 1, 2. Notice that the left-hand side of the above equation, i.e., the 
tensor H AB introduced here, depends merely on the 3-dimensional Riemannian 
metric h a b fSJ . Since in the present coordinate system the only unspecified 
component of h a t, is given by the function z(p, 9), it follows that the only unknown 
function entering into H AB is z(p,9). This property will be useful in the study 
of the integrability conditions of the field equations, since they will reduce to a 
system of equations for the single unknown function z(p, 9) [To] . 

Analogously to the vacuum case considered in [H], in studying the above field 
equations it is convenient to consider separately two cases: one in which the 
twist, D a u, and the gradient of the norm of the Killing field, D a v , are linearly 
dependent, and the other in which they are linearly independent. 
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4.2 D a v and D a uj linearly dependent 
4.2.1 The general form of the metric 

If D a v 7^ and D a uj ^ are linearly dependent, then the functions v and uj are 
functionally related, so that one can write uj — u(v). Then using equation (|4.12|) 
one obtains 

H A B = \d A Vd B V, (4.13) 
where A, B = 1,2 and we have introduced the function V by 



d A V := iryi + (duj/dv) 2 d A v. (4.14) 

From equation ()4.13|) it follows that in any local coordinate system 

H U H 22 - H 2 l2 = 0. (4.15) 

Furthermore, the integrability condition for the existence of a function V satis- 
fying ()4.13|) can be written as 

Equations (|4.15jl and ()4.16|) are direct generalizations of the conditions obtained 
for the vacuum case in [Sj. 

Recall that H A b depends merely on the 3-dimensional metric h a b, and the only 
unknown quantity entering H a b is the function z(p,9). Thus, in the coordinate 
system (x 1 = p, x 2 = 9), the equations ()4.15|) and ()4.16j) represent a system 
of two partial differential equations for the single unknown function z(p,9). It 
is straightforward to check that the choice z = z(p) and a = is a sufficient 
condition for both equations (J4.15)) and (J4.16)) to be satisfied. It can be shown, 
however, that if the pressure of the fluid satisfies P = P{p), then this condition is 
not only sufficient, but also necessary for the existence of perfect-fluid spacetimes 
with confocal inside ellipsoidal symmetry having v and to functionally related. 

Namely, suppose that the surfaces of constant pressure are given by the coor- 
dinate ellipsoids p = const, that is, P = P(p). Because of the equation of motion 
(14. 7 j) . (if p, + P 0) this implies v = v(p), which, in turn, implies u = uj(p), 
since v and uj were required to be functionally related. In addition, by calculating 
the Ricci tensor R a b for the 3-dimensional metric h a b one obtains for the single 
non- diagonal component of R a b 

Ru = ~2u(p 2 + a 2 ) sin 909' 
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where the function u = u(p,9) is defined by equation (|4.1(J|) . With this result, 
setting a = 1, b = 2 in the field equation (|4.3|) yields du/39 = , i.e., u = u{p). 
Now using the metric corresponding to the examined case of spacetimes, one 
obtains that (j4.15J) reduces to 



a 3 



du . 
P^-2u(l-u) 



cos 2 9 + a 



du 

P^ + 2«(l-u) 



p 2 = 0. (4.18) 



Since p and 9 are independent coordinates, the term in the above equation with 
cos 9, as well as the one without cos# should vanish. Excluding the solution 
u — 1, which recovers Minkowski spacetime, and the trivial solution u — 0, one 
obtains that for the perfect-fluid spacetimes under consideration the necessary 
and sufficient condition for equation ()4.18|) to be satisfied is a = 0. With this 
condition equation ()4.16|) is also identically satisfied. 

On the other hand, a = and u = u(p) imply u = z = z(p), so the system of 
equations (J4.9j) for recovering the metric function A from the quantities used in 
the projection formalism now reduces to the single equation 

I ' M -v- 2 p 2 z- l ' 2d ^=:2R, (4.19) 



sin^ d9 dp 

where the constant R could be introduced because the left-hand side of ()4.19|) 
depends merely on the variable 9, while the right-hand side depends merely on p, 
so both sides should be equal to a constant. Then it follows that A = 2Rcos9. 
Therefore, using (|4.8jl . the line element of all stationary, axially symmetric, rigidly 
rotating perfect-fluid configurations with confocal inside ellipsoidal symmetry for 
which the gradient of the norm of the timelike Killing vector field and its twist 
are linearly dependent can be given as 

ds 2 = v{p){dt + 2Rcos9d(j)f --^-(z{p)dp 2 + p 2 {d9 2 + sin 2 ## 2 )Y (4.20) 

v{p) \ ) 

provided that the surfaces of constant pressure are given by the p = const ellip- 
soids and provided that p+P ^ [T7j. Equivalently, by performing the coordinate 
transformation —v(p)~ 1 p 2 — > p 2 and introducing the function f(p) := —z(p)/v(p) 
one may write the metric also in the form 

ds 2 = v{p){dt + 2.RCOS 9d<pf + f{p)dp 2 + p 2 {d9 2 + sin 2 9d<f) 2 ), (4.21) 

which corresponds to a special case of spacetimes with co-moving ellipsoidal sym- 
metry in the sense of Definition 3. 

The algebraic type of the metrics (J4.2(J|) and (J4.21)) is Petrov type D. Unfortu- 
nately, they are not suitable for describing the interior of compact rotating bodies 
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with a regular axis of rotation. Namely, the existence of a regular rotation axis 
would require — > as 9 — > 0, for arbitrary p. This, however, is not possible 
with the metrics (|4.20j) or (|4.21|) . therefore they cannot have a regular rotation 
axis. It turns out, instead, that they are locally rotationally symmetric and they 
belong to Class I in the Stewart-Ellis classification [18J of locally rotationally 
symmetric spacetimes. 

Next we show that these perfect-fluid spacetimes can be matched to an ex- 
ternal Taub-NUT solution, provided that there exists a P = surface. The 
line-element of the Taub-NUT metric is 

ds 2 = -F(r) (dr + 21 cos 6d(j)) 2 + F~ l (r)dr 2 + (r 2 + I 2 ) (d6 2 + sin 2 6d(j) 2 ) , (4.22) 

r 2 — 2mr — I 2 . 
F (r) = r2 + p , (4-23) 

where / and m are constants [TT]. For the purpose of the matching we assume 
that there exists a p s for which P{p s ) = and that on the p = p s surface the 
relationships r = rjt (rj = const) and r = r(p) hold. We join the two spacetimes 
by requiring the continuity of the metric components at the junction surface. In 
addition, the components of the extrinsic curvature tensor of the matching surface 
are also expected to be equal on the two sides of the surface. Then one can 
express the parameters m and I of the Taub-NUT metric, as well as the junction 
parameter rj, in terms of the functions of the interior spacetime evaluated on the 
surface P = 0. That is, there exists a matching of the examined perfect-fluid 
spacetimes to the Taub-NUT metric with 

rf = J 1 "; 18 2 (4-24) 

R 2 v s f s - pi 



m 



2Jpl - tfR? 



I = rjR (4.25) 
(p 2 - 2^ 2 + ^) (4.26) 



where v s = v(p 8 ) and f s = f(p s )- Note that matchings of locally rotationally 
symmetric spacetimes to the Taub-NUT metric were considered in [19 a as well. 



4.2.2 Examples of exact solutions 

An interesting property of the metric ()4.20j) is that if the function v(p) is known, 
then z(p) can be obtained from equation (j4.12|) by performing two quadratures. 
Then the pressure and density of the fluid are obtained from equations (|4.7jl and 
1)4.11)1 . respectively, without performing any integrations. Therefore, generating 
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new exact solutions of this type is, in fact, a simple exercise. Unfortunately, this 
way there is no guarantee that the resulting spacetime will possess a meaning- 
ful equation of state. Thus, instead of giving a bunch of unphysical new exact 
solutions, in the following we rather list some interesting examples of previously 
known perfect-fluid spacetimes which can be shown to belong to the class given 
by equations (J4.2(Jj) or (J4.21)) . We wish to emphasize, however, that the list 
of previously known perfect-fluid solutions possessing confocal inside ellipsoidal 
symmetry presented here is by no means meant to be complete, there might as 
well exist other solutions in the literature, also belonging to the examined class. 

In |2U[ l2"T] a spacetime originally given by Ferwagner [22] an d later rediscovered 
by Marklund [22] is extensively studied. This spacetime represents an incompress- 
ible, rigidly rotating, Petrov type D perfect-fluid configuration. By performing 
the coordinate transformation Rsinx — ► p, dt — > dt + (R 2 / p 2 ){l — p 2 / R 2 )~^dp 
on the metric given in [20j . one arrives at 

o 4 1 

ds 2 = —^-(dt + 2i?cos 9d<p) 2 + TTT^dp 2 + p 2 (d9 2 + sin 2 9d4> 2 ), (4.27) 

it 4 1 — p 2 /R 2 

where R is a constant. The above line element has exactly the same form as 
fl4.21j) . with v(p) = p 4 /R 4 and f(p) = (1 — p 2 /R 2 )~ 1 . The density and pressure 
of the fluid are given by 




(4.28) 



Next, it has been shown in [21] that if u = uj(v) and p + 3P ^ 0, there exists 
such a coordinate system (t,x x ,x 2 in which the line element of all stationary, 
axisymmetric, rigidly rotating perfect-fluid spacetimes can be written as 

ds 2 = vix 1 ) (dt 2 + A(x 2 )d<j) 2 ) + - JL^n (x 1 ) (dx 1 ) 2 + h 22 (x 1 )[(dx 2 ) 2 + e 2 (x 2 )d<f) 2 }}, 

(4.29) 

where all possible forms of the functions e(x 2 ) and A(x 2 ) are listed in [21]. By 
setting x 1 = p, x 2 = 9, h 22 = p 2 and with the allowed choice e 2 = sin 2 9, 
which implies A = 2/2 cos one recovers the general form (|4.2U|) of the studied 
ellipsoidal spacetimes. An exact solution of the examined type was also given 
in [21], which, because of its complexity, we do not list here but rather refer 
the reader to equations (49)-(53) in the original article [21]. Finally, in [23] by 
using the 3-dimensional spin coefficient method the authors have arrived to a line 
element which has, up to a coordinate transformation, the same general form 
as that given by equation (J4.2(Jj) . but they could obtain an exact non- vacuum 
solution only for the dust case. 
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Note that for the special case of rigidly rotating dust one has P = and 
v = const. Assuming that the metric has the general form given by (|4.2U|) . it 
turns out that the only rigidly rotating dust solution satisfying the conditions 
(I4.15|) and (j4.16|) is given by z(p) = p 2 /(p 2 + v 2 R 2 ). This solution is identical 
(up to a coordinate transformation and up to renaming the constants) to a dust 
solution given in [2*5] . 



4.3 D a v and D a u linearly independent — an exact solution 

When D a v and D a u are linearly independent, the general form of the metric is 
given by (|4.8|) . In this case the components of the tensor H AB , defined by (|4.12|) . 
may be considered as components of a non-singular Riemannian 2-metric. It 
has been demonstrated in ^3] that the necessary and sufficient condition for the 
existence of the functions v (p, 9) and uj(p, 9) satisfying (j4.12j) is that the Gaussian 
curvature of this 2-metric is exactly minus one. In the present case this condition 
represents a non-linear partial differential equation for the function z(p,9). It 
can be shown that the particular choice 



p 2 + a 2 cos 2 9 p 2 
z (p, e ) = YT~2 2 2> ( 43 °) 

p z + GT 77T — p z 



where a and m are some constants, satisfies the mentioned integrability condition. 
With this choice for z, the metric ()4.8|) satisfies the full set of field equations only 
if the remaining metric functions are given by 

n Sill ^ f) 

v(p, 9) = -C\p 2 + a 2 cos 2 9), A(p, 9) = ^ + ^ — - y (4.31) 

where C is a constant. With these functions the metric of spacetime acquires the 
remarkably simple form 

ds 2 = -C 2 {p 2 + a 2 cos 2 9)dt 2 + 2a sin 2 9dtd<p + 
1 / p 2 dp 2 



C 2 \(p 2 + a 2 )(m 2 - p 2 ) 



+ d9 2 + sin 9 2 d(f) 2 ^ . (4.32) 



It is interesting to note that when the Kerr metric is recast into the form ()4.8|) 
(see |1]) then the function z(p, 9) for the Kerr metric is given by 

= ' + f"f ' / , (4-33) 
p z + a A pr — m z 

which differs from (14.30)1 merely by a negative sign factor. 
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Expression (|4.32j) represents an exact solution of Einstein's equations corre- 
sponding to a stationary, axially symmetric, perfect-fluid spacetime with confocal 
inside ellipsoidal symmetry, which does not belong to the class of solutions given 
in the previous subsection. For 9 —>■ one has —>■ which is compatible with 
the existence of a regular rotation axis. From equation (j4.11|) one obtains the 
pressure of the fluid, P. Then, using the equation of motion (|4.7|) . the density of 
the fluid, /i, can also be determined. For the particular metric (|4.32j) it turns out 
that p = — P = C 2 , i.e., the solution can be re-interpreted as a vacuum solution 
with a non-zero cosmological constant A = C 2 . We emphasize, however, that 
other perfect-fluid metrics (i.e., different from ()4.32j0 . satisfying the integrability 
conditions for spacetimes with D a v and D a uj linearly independent and possessing 
confocal inside ellipsoidal symmetry might exist as well. Proving the existence 
(or non-existence) of such metrics is still an open issue. 

5 Conclusion 

In this paper a framework has been presented which can be useful in the general- 
relativistic study of ellipsoidal shapes. This framework relies on the notion of 
a Riemannian 3-space that can be filled in with a one-parameter congruence of 
concentric ellipsoids, by the use of which one can define a conformally ellipsoidal 
spacetime. Since in general relativity the description of any shape depends on 
the relative motion of the observers performing the description, in order to carry 
out a general-relativistic study of ellipsoidal shapes one first has to specify a 
choice of preferred observers. Different choices give rise to the definitions of, 
e.g., spacetimes with co- moving ellipsoidal symmetry and spacetimes with inside 
ellipsoidal symmetry. 

An attempt to apply the presented framework in the study of compact, rapidly 
rotating bodies has been made. In particular, the applicability of the framework 
was tested by a detailed study of stationary, axially symmetric, rigidly rotating 
perfect-fluid spacetimes with confocal inside ellipsoidal symmetry. In the subcase 
when the twist, D a u, and the gradient of the norm of the time-like Killing vector 
field, D a v , are linearly related the angular dependence of the metric has been de- 
termined under the assumption that the surfaces of constant pressure are given 
by the p = const "ellipsoids" . The obtained general form of the metric (see equa- 
tions (j4.21|) and ()4.20j0 contains two unspecified functions of the coordinate p and 
it represents a Petrov type D, locally rotationally symmetric fluid configuration. 
Since this metric cannot have a regular axis of rotation, it should not be expected 
to describe the gravitational field of a compact rotating body, but rather, as it 
has been shown, interior NUT spacetimes. Some previously known exact solu- 
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tions belonging to this class have been recalled, and it has been indicated that 
new solutions of this type can be easily constructed, although their physical sig- 
nificance is questionable. In the complementary subcase in which D a u and D a v 
are linearly independent a possibly new exact perfect-fluid solution possessing a 
regular axis of rotation has been obtained. Unfortunately, it has the equation of 
state fj, = —P = C 2 and therefore corresponds to a vacuum configuration with a 
non-zero cosmo logical constant. 

The question whether in the latter subcase one might obtain perfect-fluid so- 
lutions appropriate for describing the interior of rapidly rotating compact objects 
with a regular rotation axis and a physically meaningful equation of state is still 
an open issue. Furthermore, it would be instructive to check whether the results 
obtained here still remain valid if one drops the assumption that the surfaces 
of constant pressure are given by the p = const surfaces and assumes merely 
that the surface P = is given by a coordinate ellipsoid. An additional step 
forward along the line of research presented here would be to alter the choice of 
preferred observers and to consider, e.g., spacetimes with co-moving instead of 
inside ellipsoidal symmetry. These issues are the subject of ongoing work. 
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Appendix 

A Ellipsoids in flat space 

This appendix contains a review of the properties of ellipsoids in flat Euclidean 
space on which Definition 1 in the main part of this paper is based. Note that the 
definition given in also relies on the properties of Euclidean ellipsoids listed in 
this Appendix. 

In jS] an analysis of a 3-dimensional Euclidean space filled with a congruence 
of ellipsoids of revolution with a common center and common axis of symmetry 
was given. In order to remain inside the class of axially symmetric configurations, 
the following problem was left out from |3, and left for further considerations: 
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Problem I. Make an analysis for nonsymmetric ellipsoids, analogous to that 
given for symmetric ellipsoids in 

To consider this problem, imagine a 3-dimensional Euclidean space, R 3 , filled 
with a congruence of ellipsoids with a common center. Let one of the semiaxes 
of the ellipsoids be taken as the coordinate p in space, so that p = const on a 
fixed ellipsoid. Let <fi be the angle measured around this axis, and let 9 be a third 
coordinate defined at will. Analogously to the symmetric case considered in jS], 
we define (p, 9, <fi) in terms of Cartesian coordinates (x, y, z) as follows: 



x = x (p) siia9 cos< 
y = y (p) sin6>sin</ 
z = p cos 9, 



(A.1) 



where xo(p) and yo(p) are arbitrary functions whose values equal the lengths 
of the other two semiaxes of the p = const ellipsoid. In these coordinates one 
obtains the line element ds 2 E := dx 2 + dy 2 + dz 2 of Euclidean space R 3 as: 



dst 



+ 



((x' Y cos 2 <p + (y' o y sin 2 <f>) sin 2 6 + cos 2 9 



■dp 2 



[xl cos 2 cf) + yl sin 2 <p) cos 2 9 + p 2 sin 2 9 



■dd 1 



i / 2 • 2 i , 2 2 

+ [x sin (p + y cos 



sin 2 ■ d<p 2 



+ 2{x,qx'qCos <p + yoy' sm <fi — p) cos 6 sin 9 ■ dpdO 
+ 2{—xqx' q + yoy' ) cos 4> sin sin 2 9 ■ dpd<p 
+ 2(yl — Xq) cos </>sin (f> cos9 sin # 



(A.2) 



where the prime denotes derivation with respect to p. 

The special case in which the congruence consists of confocal (oblate) ellipsoids 
is recovered by setting x\ = p 2 + a 2 , y 2 = p 2 + b 2 , where a and b are constants. On 
the other hand, for the special case of rotationally invariant ellipsoids, we have 
Xo = Ho ='■ 9, which reproduces equation (3.2) of 0: 



ds 2 E 



= ((g) 2 sin 2 9 + cos 2 9)dp 2 + 2(gg' - p) sin 9 cos 9dpd9 

i / 2 2 n i 2 • 2 n\ jn2 , 2 • 2/iji2 
+ (g COS 9 -\- p Sill 1 ' 



2 , 2-2 

+ g sin 



(A.3) 



Finally, if one considers a congruence of rotationally invariant confocal oblate 
ellipsoids, one has g 2 = p 2 + a 2 and the metric becomes 



2 p 2 + a 2 COS 2 9 2 2 2 2 

= ^ — ^ — «P + (p + a cos 

p z + a z 



+ (p 2 + a 2 ) sin 5 



(A.4) 
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Note that by letting p = const in any of the above metrics one obtains the 
metric induced on a single ellipsoidal surface with the "radius" p. By using 
the coordinate system of (jA.2|) it is straightforward to show that the extrinsic 
curvature tensor K a }> of any one of the p = const surfaces is conformal to the 
metric of the unit sphere, jab, which may be given as ds 2 {^ a b) = d6 2 + sinOdcf) 2 . 
That is 

Kab = kj ab , (A. 5) 

where k is a conformal factor, which is a function of the coordinates chosen 
to represent 3-dimensional Euclidean space. On the other hand, the intrinsic 
geometry of a singe surface does not depend on the properties of the congruence 
in which it is embedded in. Namely, in Euclidean space the surface of any p = 
const ellipsoid is conformal to the unit sphere, which implies that one can always 
introduce such coordinates on the surface in which the metric of the surface can 
be put into the form K^ a b, where K is again some conformal factor. This is an 
intrinsic property of each p = const ellipsoid. In this respect, equation (|A.5|) 
provides a connection between the "intrinsic" and "extrinsic" properties of the 
ellipsoidal surfaces in Euclidean space. 
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